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Abstract. We discuss the structure and main features of the nonlinear evolution equation proposed by this author as the 
fundamental dynamical law within the framework of Quantum Thermodynamics. The nonlinear equation generates a dynam- 
ical group providing a unique deterministic description of irreversible, conservative relaxation towards equilibrium from any 
non-equilibrium state, and satisfies a very restrictive stability requirement equivalent to Hatsopoulos-Keenan statement of the 
second law of thermodynamics. Here, we emphasize its mathematical structure and its applicability also within other contexts, 
such as Classical and Quantum Statistical Mechanics, and Information Theory. 

Keywords: Nonlinear Dynamics Dynamics; Irreversibility; Entropy; Information Theory; Quantum Thermodynamics 
PACS: 03.65.Ta. ll.lO.Lm, 04.60.-m, 05.45.-a 

INTRODUCTION 

The problem of understanding entropy and irreversibility has been tackled by a large number of preeminent scientists 
during the past century. Schools of thought have formed and flourished around different perspectives of the problem. 
But a definitive solution has yet to be found. 

We address a mathematical problem very relevant to the question of nonequilibrium and irreversibility, namely, 
that of "designing" a general evolution equation capable of describing irreversible but conservative relaxation towards 
equilibrium. Our objective is to present an interesting mathematical solution to this "design" problem, namely, a new 
nonlinear evolution equation that satisfies a set of very stringent, relevant requirements 1 1]. 

We first define four essentially different contexts within which the new equation is clearly relevant, with entirely 
different interpretations: Classical Statistical Mechanics (CSM), Classical Information Theory (CIT), Quantum Statis- 
tical Mechanics (QSM), Quantum Thermodynamics (QT). Then, we list the "design specifications" that we intend to 
impose on the desired evolution equation. We review some useful well-known mathematics involving Gram determi- 
nants and, finally, present our nonlinear evolution equation which meets the stringent design specifications. 

The discussion here is purposely devoided of our school's quite unorthodox perspective, views and hypotheses on 
the physical meaning of entropy and irreversibility. This is not only due to obvious space limitations, but mainly 
because we feel that the proposed nonlinear equation constitutes an important mathematical-physics advance in itself 
not only in the Quantum Thermodynamics contexts for which it was originally designed and developed ij^-IHHIllSlitl^ 
but also in other contexts, such as CSM, CIT, QSM, as well as possibly in Quantum Information, Biology, Sociology 
and Economics. 



FRAMEWORK A: CLASSICAL STATISTICAL MECHANICS 

Let n be a phase space, and ^ the set of real, square-integrable functions A, B, ... on £1, equipped with the inner 
product {■]■) defined by 

{A\B) = Tr(AB) = j^AB dQ. (1) 

where Tr( ) in this framework denotes Jo^ dQ.. We denote by ^ the subset of all nonnegative-definite, normalized 
functions (distributions) p in i.e., 

^ = {p in J^IP >0,Tr(p) =^pt/n= 1} (2) 

We will then consider a set {H,Ni ,.. . ,Nr} of functions in 



In Classical Statistical Mechanics, p is the Gibbs density-of-phase distribution function which represents the index 
of statistics from a generally heterogeneous ensemble of identical systems (with associated phase space fi) distributed 
over a range of possible classical mechanical states (the support of p). H is the Hamiltonian function, and A^, the 
number-of -particle function for particles of type /. 

FRAMEWORK B: CLASSICAL INFORMATION THEORY 

Let be the set of all nxn real, diagonal matrixes A = diag(flj), B = diag{bj), . . . ( n < °° ), equipped with the 
inner product (-I ) defined by 

iA\B)=TT{AB)^rj=iajbj (3) 
We denote by ^ the subset of all nonnegative-definite, unit-trace matrixes p in i.e., 

^ = {p = diag(;,,) I > 0, Tr(p) = I^LiP; = 1 } (4) 

Later we consider a set {H,Ni ,Nr} of diagonal matrixes H = diag(ey), A'l = diag(niy), . . . , N,- = diag(nr;) in Jff. 

In Information Theory fcTJ, p = diag{pj) represents the probability assignment to a set of n events, pj being the 
probability of occurrence of the j-th event. H, Ni, . . . , Nr sue characteristic features of the events in the set, taking on 
the values ej, n\j, . . . , iirj, respectively, for the j-th event. 

FRAMEWORK C: QUANTUM STATISTICAL MECHANICS 

Let be a Hilbert space (dim < °°), and Jiff the set of all linear operators A, B, ... on Jf, equipped with the 
real inner product (• | •) defined by 

(A|B)=Tr(A'^B + BU)/2 (5) 

where denotes the adjoint of operator A and Tr( ) the trace functional. We denote by the set of all self-adjoint, 
nonnegative-definite, unit-trace operators p in i.e., 

^ = {pinJ^|p'^=p,p>0,Trp = l} (6) 

We will then consider a set {H,Ni ,Nr} of self-adjoint operators in Jff, where each Ni commutes with H, i.e., is 
such that HNi = NjH, for / = 1 , . . . , r. 

In Quantum Statistical Mechanics, p is the von Neumann statistical or density operator which represents the index 
of statistics from a generally heterogeneous ensemble of identical systems (same Hilbert space J^f and operators 
{H,Ni ,Nr}) distributed over a range of generally different quantum mechanical states. If each individual member 
of the ensemble is isolated and uncorrected from the rest of the universe, its state is described according to Quantum 
Mechanics by an idempotent density operator (p^ = p = = ^-^^), i.e., a projection operator onto the span of 
some vector | y/) in J^. If the ensemble is heterogeneous, its individual member systems may be in different states, 
P\Vi)^P\Hi2)'^^^ SO on. 

H is the Hamiltonian operator, and operator A^,, for / = 1, . . . , r, is the number operator for particles of type / in the 
system (if the system has a fixed number n, of particles of type /, then Ni — n,7, where I is the identity operator on Jff). 

FRAMEWORK D: QUANTUM THERMODYNAMICS 

In our formulation of Quantum Thermodynamics HI B HI |^, the mathematical framework is the same as that 
just summarized for QSM, but the fundamental difference is in the physical meaning of the density operator. Indeed, 
QT assumes that the true individual quantum state of a system isolated and uncorrected from the rest of the universe 
is represented by a density operator p which are not necessarily idempotent. Over the set of idempotent p's, QT 
coincides with Quantum Mechanics, but it differs fundamentally from it because it assumes a broader set of possible 
states, corresponding to the set of non-idempotent p's. This way, the functional S{p) (defined in the next section) 
describes in QT an intrinsic (non-statistical) fundamental state property of the individual system. This is very different 
from the meaning that the von Neumann functional S{p) has in QSM, where it measures the degree of heterogeneity 
of the ensemble whose statistics are represented by p. 



MEAN VALUE FUNCTIONALS AND 5-FUNCTIONAL 



From here on, our notation allows us to treat at once the four contexts just defined. For reasons to become apparent 
below, the elements H, Ni,. . . , Nr introduced in either context, will be called the generators of the motion (MG). We 
assume that such sets always contain at least element H, that we call the Hamiltonian MG. 

For each MG, we then define a mean value functional on ^ as follows 



Moreover, we define the 5-functional |(^ on ,^3^ as 



tn{p;Ni) = TrpNi = (^/p|v^A',■) 



Sip) = -m(plnp) = -ki^l^lnp) 



(7) 



(8) 



Depending on the context, the 5-functional represents the statistical uncertainty as to the actual classical or quantum 
state of a system, the information carried by the occurrence of one of the possible events (or the degree of uncertainty 
as to which will be the next event), or the thermodynamic entropy. 

For each given set of values (//), (A^i), . . . , (A^,-), in the range of the mean value functionals (Eqs.|7]l corresponding 
to the GM's, we consider the subset of all elements p in ^ that share the given mean values, i.e., 

^m,{N,):-,{Nr)} = iP -^I'HP^H) = {H),m{p-M) = {Ni) for / = 1, . . . ,r} (9) 

On each such subset, i.e., for fixed mean values (//), (A^i), . . . , {Nr) of the generators of the motion, the 5-functional 
(Eq.[8]i achieves a unique maximum at the point 



p = g-" exp ( - j3// + v,A^, ) where a = InTr 



(10) 



exp(^-j3// + L[^iV,M^ 
Vi {{H),{Ni),. . . , {Nr))- It is noteworthy that the maximum-5 

(11) 



and, of course, /3 = j3 ((//>, (M>, ... , {N-)) and V/ 
points satisfy the condition 

Vplnp = -a^-b^H + YJi^^Ci^Ni 

for some real numbers a, b and c,-, / = 1, . . . ,r. In words, the maximum-S' element p is such that ^/plnp lies in the 
linear manifold generated by elements ^/p, \fpH, yfpN\ , . . . , Nr- Condition[TT|is satisfied not only by the maximum-S' 
elements given by Eq.[TOl but also by the elements given by 



p=e'"BQxp{-bH + YIi=i CiNi ] B where a = In Tr 



Bexp(-bH + l,'-^^CiNi]B 



(12) 

where b = b{B;{H) ,{Ni) , - - - ,{Nr)), Ci = c,- (B; (//), (A^i ),..., (A^,-) ) and B is any idempotent element in ^ (i.e., 
B^ — B)- Clearly, Eq. [12] reduces to Eq. [TO]iff B = I (I — constant function equal to 1 on the whole Q. in CSM; 
/ = diag(l ) in CIT; / = identity operator on ^ in QSM and QT). 



DYNAMICAL LAW DESIGN SPECIFICATIONS 

Our scope is to design a function F{-) such that every solution p (?) of the autonomous differential equation 

i-p{t)^F{p{t)) (13) 

with p(0) anywhere in ^ satisfies the following conditions for all f 's, — oo < f < oo; 

(i) p(f ) lies entirely in 3^ (no forward nor backward escape times); 

(ii) m{p{t);H)= m{p{0);H), and m{p{t);Ni) = m{p{0);Ni) for / = 1, . . . ,r; 

(iii) S{p{t + u))>S{p (t)) for all m > 0; 

(iv) within each subset ^{{H),{Nt)....,{Nr)}' maximum-5 element given by Eq. [TO]i is the only equilibrium solution 
that is stable according to Lyapunov 121; all other equilibrium elements must not be stable. 

Notice that requirement (iv) is most restrictive. For example, within QSM, it rules out the von Neumann evolution 
equation {F{p) — —i{Hp — pH)/h) because all the stationary density matrices ( p such that Hp = pH) are stable 
according to Lyapunov and, in general, there are many more than a single one within each set i^[(H).{Ni),...,{Nr)}- 



SOME NECESSARY MATHEMATICAL BACKGROUND 



Given a subset of elements A, B, ... ,Z in ^ , we denote by M(A,Z?, . . . ,Z) the Gram matrix 



(A|A) 


(A|B) . 


■ m) 


(B|A) 


{B\B) . 


■ m) 


{Z\A) 


{Z\B) . 


■ (z|z) 



(14) 



where ( I ) is the real symmetric inner product defined on ^ . We denote by G{A,B, . . . ,Z) the Gram determinant 
of A,B, . . . ,Z with respect to inner product (-I ), i.e., G{A,B,.. . ,Z) = det[M(A,B,. . . ,Z)]. Matrix M(A,B, .. . ,Z) is 
nonnegative definite and G{A,B, . . . ,Z) is nonnegative. Elements A,B,...,Z are linearly independent (LI) iff their 
Gram determinant G{A,B, . . . ,Z) is nonzero and, hence, strictly positive. 

Given a subset of elements A,B, . . . ,Z in .if, we denote by L{A,B, . . . ,Z) the linear manifold spanned by all linear 
combinations with real coefficients of the elements A,B, . . . ,Z. With respect to the inner product ( j ) defined on 
we denote the projection of a given element V in ^ onto a linear manifold L by the symbol {V)i^. (V )^ is the unique 
element in L such that iiV)^\X) = {V\X) for all X in L. 

The theory of Gram determinants, very seldom used in the physics literature, offers a useful explicit way of 
writing the projection (V)^ of V onto a given linear manifold L. Let the given linear manifold be L = L{A,B, . . . ,Z), 
where elements A,B, ... ,Z need not be LI. Select any subset of LI elements Ei,E2, . . . ,E,„ spanning L, i.e., such 
that G{Ei,E2,...,E„,) > and L{Ei,E2, . . . ,E,n) = L. By the definition of (V)^, {{V)^\Ej) ^ {V\Ej) for every 
y = 1,2, . . . ,OT, and (V)^ = T4L1 vtEi, where v, are real scalars. Thus, 



i:jLiVi{Ei\E,) = {V\Ej) for ./ = 1,2, . . . ,m 

Because {Ei\Ej) — [M{E\,E2, . . . ,£'m)],j and the elements £1 ,£2, . . . j-E,,, are LI, Eqs. 
the V('s to yield 



■lJ^^,{V\Ej)[M{EuE2,...,E,n)- 



for / 



1,2,. 



and, therefore. 



{V\ = rr=iL7=iiV\Ej)mEuE2,.. . ,E„,)-\Ei 
Alternatively, Cramer's rule yields the equivalent, more elegant expression 



(15) 

] are LI and can be solved for 
,m (16) 

(17) 



G{E\,E2, . . . ,E„, 
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{Efji 1 Ef^ ) 



(18) 



Below, for any given element p in the set we will need to consider the projection of ^/plnp onto the linear 
manifold L{y^, ^H, y^Nr) where y^H, y^Nr are not necessarily LI. Using Eq. [18] 

and Definitions [1] [3] and |5] of the inner product ( j ), we find 



(VPlnP)L(yp,yp//,ypw, 



v/pA',-: 



X det 



Ti-pRl 





Trplnp 
TipRilnp ^Tip{Ri,Ro} 

Trp/?,lnp 



2' 

^Ttp{R„Ro} 



G{^Rq,^Ru---VpRz) 



^Trp{Ro,Ri} 
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(19) 



(20) 



where {A,B} = AB + BA, and Ro,Ri, . . . ,R^ are a subset of elements in ^ such that L(y/p7?o, ^/pRi , 
L( VP, VpH, VPM , . . . , VPM ) and G{^Rq, VP^i , . . . , ^R,) > 0. 



•,VP^z) = 



STEPEEST-5-ASCENT NONLINEAR EVOLUTION EQUATION 



With the above background, the nonlinear evolution equation proposed by the author H El H [lltl to meet our 
design specifications takes the compact form: 

^p(0 = -j^[H,p] - ^ {^D{p)+D\p)^) (21) 

D{p) = Vplnp - (Vplnp)^(^^^^^^ ^^_.) (22) 

where [i/,p] = Hp — pH (= within CSM and CIT), fi is the reduced Planck constant (playing a role only within 
QSM and QT), t is a characteristic time constant, H,N\ ,Nr are fixed GM's. 

Notice that D(p (f )) is orthogonal to the linear manifold spanned by ^{t),^/p{t)H , y/p{t)Ni y/p{t)Nr and the 
term i[H,p] is orthogonal to y/pD{p) + D'^ (p)^. 

Within QT, if the isolated and uncorrected system AB is composed of subsystems A (Alice) and B (Bob) that are 
either interacting (H — H^^^I^ + correlated (S ^ ^^(Xi/^ + where S — ~kB {Prwp ) In p ), or both, 

the proposed equation takes the form ||6lll2|l (we assume for simplicity that A and B have no non-Hamiltonian GM's): 



-£ = -^[H,p] + -—{^D^+dI^) ®p^ + _-p^® (^D^+4VP^) (23) 

{Hf = Tr,[{I^®p^)H] [Hf = Tr^[(p^®4)//] (24) 

{Sf = Tr^ [(4®Pb)5] [Sf = Tr^ [{Pa®Ib)S] 

Despite the nonlinearity, the structure of the non-Hamiltonian terms in the equation prevents "no-signaling" violations. 

9 dp / p , 

All zero entropy states (p = p), even if A and B are entangled, obey the Schroedinger equation — = — — [H,p\, 

at n 

thus not contradicting any of the results of QM, although within QT, these solutions, including the stationary states of 
QM (p^ = p, pH ~ Hp) are very weakly unstable limit cycles or equilibrium states 0131. 



The proofs that this equations satisfy our design specification are in Several other intriguing 

features including Onsager's reciprocal relations are discussed in Refs. IToi llU Fl ill . Eq. |23] generalizes Eq. 1211 to 



es £,q. j/. 

composite systems consistently with the additional design specs required to avoid non-locality paradoxes 1^61 [l2ll . 

We believe that in view of its intriguing properties, this evolution equation constitutes an important mathematical 
"tool" for a variety of non-equilibrium relaxation problems, not only within our QT, but also within different contexts, 
such as CSM, CIT, and QSM, as discussed here, as well as Quantum Information, Biology, Sociology and Economics. 
Even in the Thermodynamics context, where different schools of thought notoriously have contrasting perspectives on 
the physical meaning of entropy and irreversibility, many important insights can nevertheless be extracted from the 
richness of structure and the well-behaved and self-consistent features of this relevant nonlinear dynamical equation. 
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